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ABSTRACT 

After  some  discussion  of  waU  fires,  the  effects  that  a ceiling  has  on  a wall  fire  are  analyzed  and 
discussed  qualitatively.  There  are  two  kinds  of  effects:  first,  when  the  upper  sill  of  the 
ventilating  opening  lies  below  the  ceiling  (as  is  usual),  a layer  of  hot  gas  is  trapped,  which  heats 
the  waUs  by  convection  and  radiation;  the  heated  ceiling  also  radiates  to  the  walls.  A calculation 
which  uses  the  wall-fire  computer  model  SPREAD  is  carried  out  to  demonstrate  how  these 
effects  can  be  calculated  now.  Second,  when  the  flame  impinges  on  the  ceiling,  it  bifurcates  and 
spreads  horizontally,  rather  than  vertically.  A simple  algorithm  is  presented  which  calculates 
the  spread  rate  of  these  horizontal  flame  extensions,  as  well  as  the  modified  pyrolysis  rate,  due 
to  the  modified  radiation  feedbacks. 


Key  words:  algorithm;  ceilings;  ceiling  fires;  computer  models;  fires;  fire  growth;  fire  spread; 
SPREAD;  wall  fires 
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1.  INTRODUCTION 


Wall  fires  are  significant  for  two  related  reasons:  first,  because  fires  spread  upward  swiftly  on 
vertical  surfaces,  often  leading  to  rapid  flashover  in  rooms.  Second,  the  upward  spread  rate, 
especially  in  a comer,  is  a sensitive  function  of  the  fiammability  of  the  material.  Therefore  it 
can  be  used  to  estimate  its  fiammability.  Indeed,  full-scale  comer  tests  have  been  developed  to 
assess  the  fire  hazards  of  material  linings. 

Accurate  modeling  of  wall  fires  will  enable  the  flammability  and  hazard  of  wall-covering 
materials  to  be  assessed  from  bench-scale  tests  of  the  materials,  rather  than  having  to  rely  on 
more  expensive  and  time-consuming  full-scale  tests.  Also,  such  a modeling  capability  will  allow 
the  determination  of  the  interactions  among  fuel  packets  (e.g.,  items  of  furniture),  the  enclosure, 
and  flame  spread  along  combustible  walls  and  ceilings. 

Our  inability  to  predict  these  interactions  accurately  has  led  to  difficulties,  such  as  confusion  in 
the  interpretation  of  European  and  American  comer  room  fire  tests.  The  EUREFIC  test  (ISO 
975;  see  Proceedings,  1991),  which  evaluates  interior  lining  materials,  is  a comer  test  where  the 
walls  and  ceiling  are  both  lined  with  the  same  material.  The  comer  tests  used  in  the  USA 
(ASTM;  UBC  Standard  42-2)  do  not  have  the  ceiling  lined.  The  results  can  be  very  different 
with  the  two  tests  (Sumathipala,  1993).  It  would  be  very  useful  to  come  up  with  a model  which 
can  cover  both  cases. 


2.  BACKGROUND 

A number  of  studies  of  wall  fires  have  been  made,  and  several  wall-fire  models  created  to 
simulate  the  ignition  and  upward  spread  of  wall  fires.  For  studies,  see  Ahmad  and  Faeth  (1979), 
Hasemi  (1986),  Orloff  et  al  (1975,  1977,  1979),  for  example.  For  models,  see  Saito  et  al 
(1985),  Mitler  (1986,  1988,  1990),  Quintiere  et  al  (1986),  Delichatsios  (1987),  Chen  et  al 
(1993),  Brehob  and  Kulkami  (1993),  etc.  NIST  has  developed  reasonably  successful  models 
which  simulate/predict  wall  fires  and  upward  flame  spread  under  unconfined  (no-ceiling) 
conditions  (Mitler,  1989,  1990;  Steckler  and  Mitler,  1988).  The  model  name  is  SPREAD. 

Most  of  these  studies,  however,  have  ignored  the  effect(s)  of  the  ceiling.  Thus,  most  extant 
models  assume  (in  effect)  that  the  wall  is  infinitely  high.  In  fact,  the  presence  of  a ceiling  has 
several  important  consequences,  and  the  effect(s)  of  the  ceiling  on  room  fires  is  one  of  the  most 
notable  missing  pieces  in  our  ability  to  accurately  predict  the  development  of  room  fires. 

Friedman  (1975)  reported  that  radiative  feedback  from  a hood  tripled  the  burning  rate  of  a 
PMMA  pool  fire,  compared  to  its  burning  in  the  open.  Hinkley  et  al  (1969-84)  examined  flames 
spreading  under  a corridor  ceiling,  and  their  substantial  effect  on  the  ignition  of  floor  samples, 
among  other  things.  Alpert  et  al  (1981)  studied  wall  fires  spreading  into  a corridor-like  ceiling 
channel,  similar  to  Hinkley’s.  These  experiments  demonstrated  that  wall  burning  rates  in 
confined  spaces  increased  significantly  over  open-burning  rates.  While  large-scale  experiments 
on  flame  extensions  has  been  done  by  Gross  (1989)  for  axisymmetric  and  comer  fires  with 
noncombustible  ceilings,  and  by  Kokkala  (1991)  for  axisymmetric  fires,  no  work  has  been 
reported  for  extensions  against  a flat  wall  with  a ceiling. 
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Mitler  and  Steckler  (1993)  reviewed  the  literature,  quantified  the  effect  (in  part),  and 
demonstrated  the  importance  of  ^he  ceiling  on  the  development  of  enclosure  fires.  The 
conclusion  was  that  ceilings  do  sed  have  significant  effects  on  the  development  of  fires  in 
enclosures.  Thus,  as  has  been  demonstrated  by  these  authors,  the  effects  of  ceilings  on  the 
progress  of  a compartment  fire  are  substantial,  and  need  to  be  included  in  any  compartment 
waU-fire  calculations.  Although  the  available  studies  address  various  aspects  of  the  ceiling-fire 
problem,  they  collectively  do  not  provide  a minimally  complete  set  of  data  (e.g.,  flame  shape, 
thickness,  temperature  distribution,  heat  flux,  etc.),  for  the  development  and  validation  of  a wall- 
ceiling or  comer-ceiling  fire  algorithm. 


3.  EFFECTS  OF  THE  PRESENCE  OF  A CEILING 
3. 1 Flame  Tips  Below  the  Ceiling 

The  most  obvious  effect  of  the  presence  of  a ceiling  is  the  development  of  a hot,  smoky,  and 
oxygen-vitiated  layer  under  the  ceiling.  This  layer  increases  the  heat  fluxes  to  various  surfaces 
in  the  room.  The  added  heat  flux  to  the  already-pyrolyzing  regions  of  a burning  waU  accelerates 
the  burning  and  spread  rates.  The  added  heat  flux  to  the  preheating  region  accelerates  the  spread 
rate  still  further.  This  is  illustrated  by  the  schematic  drawing  in  figure  1.  This  occurs  even 
before  the  flame  "splashes"  against  the  ceiling;  when  the  flame  is  longer,  additional  effects  come 
into  play,  as  described  in  the  next  subsection. 


3.2  Flame  Length  Greater  than  the  Wall  Height 

When  the  wall  flame  reaches  the  ceiling,  it  splashes,  and  the  flame  "extension"  is  forced  to 
move  at  right  angles  to  its  previous  trajectory.  Part  of  the  flame  extension  may  be  an  oval  or 
other  shape,  projected  on  the  ceiling,  as  shown  in  figures  2 and  3.  Because  of  that,  the  view 
factors  for  radiation  feedback  from  this  part  of  the  flame  to  both  the  pyrolyzing  and  the 
preheating  sections  of  the  wall  are  increased;  this  is  indicated  in  figure  4.  This  enhancing  effect 
is  in  addition  to  the  extra  flux  from  the  hot  gas  layer,  discussed  above.  The  burning  as  well  as 
the  upward  spread  rates  of  the  wall  fire  will  thus  be  further  accelerated. 

As  soon  as  the  pyrolysis  front  reaches  the  ceiling,  "upward"  spread  must  be  reinterpreted  as 
concurrent-flow  spread  (in  the  horizontal  direction),  and  the  rates  will  depend  on  the  geometry 
and  configuration  of  the  enclosure  and  the  resulting  flows  and  embedded  flames.  These  can  be 
quantitatively  calculated,  once  we  know  the  length  and  shape  of  the  flame  extensions. 

Experiments  have  demonstrated  that  the  flame  moves  out  along  the  wall/ceiling  intersection,  as 
shown  in  figure  3,  rather  than  continuing  to  move  in  a direction  as  close  as  possible  to  the 
previous  (upward)  direction,  as  shown  in  figure  2,  or  fanning  out  across  the  ceiling  as  a 
semicircle,  as  is  indicated  in  figure  5.  The  reason  it  does  this  is  probably  because  of  the 
formation  of  vortices,  as  first  hypothesized  by  Williamson  et  al  (1984),  and  demonstrated  for 
the  comer  fire  by  Satoh  (1989). 
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Quintiere  (1988)  described  an  experiment  carried  out  by  B.  Lee  in  1988  (see  Figure  6).  The 
back  wall  of  an  ASTM  enclosure  (2.44x3.56x2.44  m --  i.e.,  8x12x8  feet  --  see  Anon.,  1983, 
1993)  was  partially  covered  by  6.35  mm-thick  (1/4")  sheets  of  PMMA  in  the  shape  of  a "T." 
The  bottom  was  1.22  m (4’)  wide,  and  the  cross-bar  above  was  1.22x2. 44m  (4’x8’);  (the  results 
would  probably  not  have  been  any  different  if  the  entire  wall  had  been  covered).  It  is  important 
to  note  the  flaming  "arms"  which  develop  along  the  wall-ceiling  intersection  (they  appear  to  lie 
principally  against  the  wall).  These  arms,  or  "extensions,"  have  the  length  and  width  d (the 
thickness  will  be  discussed  later)  ~ see  Figure  7.  Predicting  the  further  rate  of  spread  of  the 
fire  depends  in  part  on  knowing  d and  £^.  In  order  to  calculate  the  length  and  shape  of  the  flame 
extensions,  the  new  view  factors,  and  the  consequent  accelerated  preheating  and  subsequent 
advance  of  the  pyrolysis  front,  it  is  necessary  to  understand  quantitatively  the  flow  of  the  thermal 
plume. 

Cleary  and  Quintiere  (1991)  have  developed  a wall/ceiling-fire  model  based  on  the  Lee 
experiments.  They  simplify  the  problem  by  assuming  that  the  width  5 of  the  arms  is  w/2,  where 
w is  the  width  of  the  wall  fire;  they  also  assume  that  there  is,  additionally,  a semicircular  flame 
zone  against  the  ceiling,  as  shown  in  figure  5.  This  formulation  is  limited,  and  leads  to  at  least 
one  quantitatively  incorrect  conclusion:  if  w > 2h^,  then  5 > h^.  That  is,  the  width  of  the 
flame  extension  is  greater  than  the  height  of  the  wall. 


3.3  Ceiling  Ignition  and  Downward  Flame  Spread 

In  order  to  model  wall/ceiling  flame  spread  more  completely,  consideration  must  also  be  given 
to  two  other  phenomena:  first,  the  effects  of  the  ignition  of  a combustible  ceiling.  Second,  the 
downward  spread  of  the  fire;  although  this  is  counter-flow  spread,  just  like  the  lateral  spread  of 
the  wall  fire,  it  has  been  observed  to  be  much  more  rapid  than  the  latter  (Karlsson  1989,  1992). 

Small-scale  experimental  work  on  burning  and  spread  rates  across  (channeled,  two-dimensional) 
combustible  ceilings  has  been  done  by  Atreya  (Atreya  and  Mekki,  1992)  for  laminar  flames  and 
by  Pello  for  turbulent  as  well  as  laminar  flames;  the  latter,  for  combustible  ceilings.  One  of 
PeUo’s  observations  is  that  the  onset  of  turbulence  shortens  the  flames,  presumably  because 
mixing  with  ambient  air  is  enhanced,  permitting  faster  and  more  efficient  combustion  (Zhou, 
1992). 

In  order  to  establish  predictive  algorithms  to  describe  the  above  effects  correctly  in  a model, 
we  need  to  develop  a quantitative  understanding  of  fluid  flows,  entrainment,  and  combustion 
processes  at  and  under  the  ceiling,  upon  the  impingement  of  the  flames  on  the  ceiling. 


3.4  Other  Effects 

Another  effect  is  one  observed  by  Saito  (1993)  for  corner  fires:  he  carried  out  a number  of 
comer  fire  experiments  with  a burner  and  inert  walls;  for  the  same  rate  of  heat  release,  the  heat 
fluxes  to  the  wall  and  the  resultant  wall  temperatures,  in  the  steady  state,  were  smaller  when  a 
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ceiling  was  placed  over  the  comer,  contrary  to  what  one’s  intuition  would  have  predicted.  His 
tentative  explanation  for  the  phenomenon  is  that  the  "back  pressure"  from  the  ceiling  slows  the 
upward  flow  of  the  thermal  plume,  resulting  in  a smaller  heat  transfer  coefficient  and  therefore 
in  a lower  convective  heat  transfer  to  the  wall,  at  least  in  the  lower  portions  of  the  flame. 


4.  THEORETICAL  CONSIDERATIONS 

4. 1 Effects  of  Trapped  Gas  Layer  and  Hot  Ceiling 

We  now  discuss  the  first  of  the  effects  discussed  in  Section  3.  This  part  of  the  ceiling  effect  can 
be  theoretically  calculated  (i.e.,  modeled)  now,  in  fact,  for  the  "early"  stages  of  a wall  fire, 
before  the  flame  has  impinged  on  the  ceiling.  The  procedure  is  straightforward,  though  (at 
present)  awkward.  This  is  done  as  follows:  the  upward  spread  of  a fire  on  a wall  of  height  h^ 
is  calculated;  it  is  done  here  with  the  program  SPREAD.  This  yields,  among  other  things,  the 
rate  of  heat  release  from  the  slab/wall.  This  power  output  is  then  used  as  an  input  to  a room- 
fire  program  (FIRST  (Mitler  and  Rockett,  1987)  was  used  for  this  purpose).  The  latter  then 
yields  the  ceiling  and  upper-layer  gas  temperatures  as  a function  of  time,  and  the  resulting 
convective  and  radiative  heat  fluxes  back  to  the  wall.  The  original  upward- spread  program 
(SPREAD,  again)  is  then  rerun  with  those  fluxes  as  input.  This  results  in  enhanced  growth  and 
spread  rates.  Note:  even  if  SPREAD  does  not  give  exact  answers,  it  will  act  as  a baseline,  and 
the  work  to  be  given  will  show,  with  reasonable  accuracy,  the  departures  from  those  baseline 
results. 

As  an  example,  curve  NC  in  figure  8 shows  the  power  output  from  a 2-foot-wide,  8-foot-high 
slab  of  PMMA,  ignited  at  the  bottom,  as  a function  of  time  (this  includes  the  RHR  from  the 
igniter).  Ignition  occurs  at  t = 198  s;  the  top  of  the  slab  is  reached  at  t = 677  s,  and  the  RHR 
there  is  315  kW.  The  RHR  continues  to  increase  slowly  thereafter.  In  curve  C,  the  results  of 
the  first  iteration  of  the  procedure  described  above  is  shown;  it  is  seen  that  ignition  occurs  about 
20  seconds  earlier,  at  t = 178,  and  that  the  upward  spread  rate  is  higher:  the  top  of  the  slab  is 
reached  about  170  seconds  earlier,  at  t = 504,  and  the  RHR  at  that  moment  is  350  kW. 
Incidentally,  that,  in  turn,  means  that  the  RHR  used  in  FIRST  should  be  the  enhanced  rates. 
Proceeding  in  this  way,  one  (should)  iterate  to  convergence. 

Note  that  a room-fire  model  uses  the  correct  ceiling  height,  whereas  SPREAD  implicitly  assumes 
a ceiling  so  high  that  the  flames  never  impinge  on  it.  The  procedure  indicated  above  will  be 
carried  out  automatically,  without  the  need  to  iterate  back  and  forth,  once  SPREAD  is  embedded 
into  a room-fire  model.  SPREAD  is  the  "test-bed"  prototype  program. 


4.2  Flow  Patterns  of  the  Thermal  Plume 

Additional  effects  arise  from  radiative  feedback  from  the  flames  under  the  ceiling.  The  dashed 
curve  in  figure  9 indicates  schematically  how  the  further  (lateral)  extension  of  the  fire  along  the 


5 


arms  on  the  wall  would  accelerate  the  RHR  curve  even  more.  In  order  to  calculate  this 
feedback,  one  must  calculate  the  view  factors.  But  in  order  to  do  that,  it  is  necessary  to  know 
the  positions  and  areas  of  the  flame  extensions.  Since  the  flames  will  follow  the  thermal  plume 
flows,  the  first  task  is  to  determine  what  the  flows  are. 

Referring  to  figure  3,  we  see  that  there  must  be  hot  gas  flow  along  the  arms.  Indeed,  it  is 
appropriate  here  to  point  out  that  the  gas  flows  are  not  intrinsically  visible;  their  existence  is 
inferred  from  simple  mass  and  momentum  conservation.  Their  shape(s),  location,  and  direction 
become  evident  only  when  they  transport  a visible  component,  such  as  particles  or  luminous 
flames. 

There  probably  is  a ceiling  jet  as  well,  though  it  is  not  evident  because  it  does  not  carry  along 
much  (if  any)  flame.  This  bifurcation  of  the  flow  into  arms  and  into  a ceiling  jet  is  expected 
to  be  approximately  independent  of  the  existence  of  combustion  in  the  gas  phase  (except  insofar 
as  combustion  gives  rise  to  buoyancy). 

For  the  moment  we  will  assume  that  the  arms  are  contiguous  with  the  wall  only.  That  is, 
except  for  the  thickness  of  the  (flame)  arms,  there  is  no  part  of  the  arms  which  is  in  contact 
with  the  ceiling.  This  is,  of  course,  subject  to  experimental  verification. 

The  upward  mass  flow  rate  in  the  plume  at  any  height  x is 

m{x)  = wv{x)pg{x)B{x)  (1) 

where  w is  the  width  of  the  plume,  Pg  its  density,  6{x)  is  the  thickness  of  the  plume  at  the  height 
X,  and  v(x)  the  (mean)  velocity  across  the  plume  at  x. 


Then  at  the  ceiling,  the  upward  mass  flow  of  the  plume  is 

where  the  subscripts  "c"  indicate  that  the  various  quantities  are  evaluated  at  the  Ceiling.  A 
fraction  F of  the  flow  m^  continues  along  the  ceiling,  away  from  the  wall  (that  is,  in  the 
direction  z,  orthogonal  to  x and  y);  we  refer  to  this  fraction,  as  before,  as  the  ceiling  jet;  we 
may  also  refer  to  it  as  "the  central  jet."  The  remainder  (the  fraction  1 - F)  goes  ofi*  along  the 
wall-ceiling  intersection.  Again  as  before,  we  call  the  latter  jets  the  "arms,"  or  "extensions" 
(see  Figure  7),  of  width  5 and  length  In  Eqs.(13)  to  (15),  below,  we  take  the  width  of  the 
central  jet  to  be  Wj,  analogous  to  6 for  the  arms.  The  analogue  to  the  arm  length  is  the 
central  jet  length  "a,"  discussed  in  the  section  on  radiation  feedback. 

We  assume  that  there  is  a negligible  loss  of  energy  from  the  turn  at  the  ceiling;  hence  the 
horizontal  flow  speed  must  be  largely  unchanged  from  the  upward  speed  at  the  ceiling,  v^. 
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Assuming  no  discontinuity  in  the  thickness  of  the  flow,  the  thickness  of  the  arms  where  they  join 
the  upward-moving  plume,  is  the  same  as  6^.  Finally,  assuming  the  density  of  the  gas 
constituting  the  arms  also  to  be  unchanged,  the  mass  flow  along  the  arms  must  be 


Thus, 


rh„  = (1  -F)/n^  = 


F = 1 - 25/w 


(4) 

(5) 


It  is  also  clear  from  Eq.(4)  and  from  symmetry  that  we  can  never  have  the  sum  of  the  widths 
of  the  arms  greater  than  the  width  of  the  plume;  that  is. 


b < w/2  (6) 

We  now  obtain  an  approximate  expression  for  5 as  a function  of  w.  Since  little  experimental 
data  on  flows  exists,  we  will  infer  as  much  as  possible  from  symmetry  arguments.  First, 
consider  the  behavior  of  the  hot  gases  across  the  ceiling,  when  the  wall  and  the  Are  are  both 
extremely  wide  (that  is,  w > > h^)  — this  is  illustrated  in  flgure  10;  the  wavy  upwardly- 
pointing  arrows  at  the  top  indicate  the  flow  of  the  hot  product  gases  (the  ceiling  jet).  It  is  clear 
from  symmetry  that  away  from  the  edges  of  the  Are,  the  flow  must  be  nearly  two-dimensional, 
as  if  the  ceiling  jet  were  traveling  down  a corridor;  thus,  F approaches  unity.  Note  that  we  have 
omitted  showing  the  arms,  in  figure  10. 

We  evidently  must  always  have  5 < h^.  Then,  assuming  that  the  thickness  of  the  arms  is  the 
same  as  the  flocal)  thickness  of  the  plume  {6^  = B^,  and  that  the  speed  of  the  flow  does  not 
change  at  the  turn,  then  for  the  case  w > 2h^,  we  must  have  b < w/2;  the  "excess"  must  travel 
along  the  ceiling,  orthogonal  to  the  arms. 

We  know  from  the  experiments  carried  out  by  B.  Lee  and  others  (figure  6)  that  there  will  be 
some  flow  down  the  two  arms,  as  shown  in  figure  7;  that  is,  F < 1.  In  light  of  that,  consider 
the  opposite  extreme:  a very  narrow  fire,  as  shown  schematically  in  figure  1 1 . Assuming  that 
the  flame  always  spreads  laterally  along  the  wall-ceiling  intersection,  there  are  two  possibilities: 

a)  besides  the  flow  along  the  two  arms,  a narrow  tendril  of  gas  continues  along  the  ceiling;  this 
may  spread  out  radially  into  a fan,  or  behave  some  other  way. 

b)  there  may  be  no  such  flow,  until  some  critical  width  w^  is  achieved.  If  that  were  so,  then 
for  w slightly  greater  than  w^,  a "tendril"  of  gas  begins  to  escape. 

Finally,  for  an  intermediate-width  fire,  we  must  have  a mixed  regime,  where  there  is  defiinitely 
some  flow  along  the  ceiling,  besides  that  along  the  arms,  as  shown  in  figure  7. 

We  now  try  to  express  these  points  in  an  explicit  relationship.  When  w is  very  narrow,  it  is 
clear  that  b must  be  small.  For  larger  w,  on  the  other  hand,  the  vertical  width  b of  the  extension 
is  (most  likely)  only  weakly  dependent  on  the  horizontal  width  w of  the  original  flame.  We 
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assume  that  to  be  the  case.  Thus,  we  might  expect  a dependence  for  5(w)  such  as  is  given  by 
the  following  simple  expression  which  satisfies  the  above  observations  for  all  three  cases  (w 
narrow,  intermediate,  and  wide): 

6(w)  * ^[l  - e'/”"] 


where  A and  p are  to  be  determined.  To  this  end,  we  note  that  Alpert  (1975)  found  that  the 
initial  ceiling-jet  thickness  for  an  axisymmetric  plume  is  5 to  8%  of  the  room  height.  Motevalli 
and  Marks  (1991)  found  that  the  thermal  boundary  layer  approaches  11.2%,  asymptotically,  for 
some  small,  axisymmetric  fires.  For  a comer  fire,  Karlsson  and  Magnusson  (1990)  found  the 
thickness  to  be  5%  of  the  room  (ceiling)  height,  while  Quintiere  (1993)  took  it  to  be  8%  of  the 
room  height.  "Room  height"  may  be  interpreted  as  "fire-surface  to  ceiling"  height. 

If  5(w)  is  given  by  Eq.(7),  then  (5)  yields 

F = \ - e-P'^)  (8) 

w 

It  is  clear  that  as  w approaches  zero,  F must  go  to  some  constant  value  which  is  nonnegative  and 
no  greater  than  1 : 

0 < 2Ap  < 1. 


If  A < l/(2p),  the  result  is  shown  as  curve  A in  figure  12.  For  the  case  that  F approaches  zero 
as  w approaches  zero,  A = l/(2p),  and 


F = 1 - 

This  is  shown  as  curve  B in  figure  12. 


pw 


(9) 


On  the  other  hand,  it  is  possible  (as  discussed  earlier)  that  until  some  critical  width  w^  is 
achieved,  there  is  no  contribution  to  the  central  jet.  If  this  is  in  fact  the  case,  i.e.,  if  F = 0 for 
w < w^j,  then  we  must  have 


5 = w/2  for  w < Wjj. 

Finally,  for  w > w^.,  we  must  have 

5 < w/2. 

That  is,  if  there  indeed  exists  a critical  width  w^,,  then 


w/2 

w < w^ 

(10) 

c 

< w/2 

W > W- 

The  simplest  such  approximation  is 
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w/2 


(11) 


d = 


We/2 


w < 
w > 


although  it  is  certainly  possible  to  have 

W(j/2  < 6 < w/2. 

This  dependence  is  shown  as  the  solid  curve  in  figure  13,  and  from  Eq.(5),  the  dependence  of 
F on  w would  be 


F = 


0 

w -w^ 
w 


w < 


w > 


(12) 


as  shown  in  curve  C,  figure  12.  The  initial  slope  of  the  solid  curve  in  Figure  13  is  0.5.  The 
dashed  curves  in  that  figure  correspond  to  Eq.(7):  the  upper  dashed  curve  corresponds  to 
Ap  > 0.5  while  the  lower  dashed  curve  corresponds  to  Ap  < 0.5. 


Since  a fraction  of  the  flow  has  been  detoured  into  the  arms,  the  initial  width  of  the  ceiling  jet 
may  be  smaller  than  w.  It  will  surely  be  so,  if  its  thickness  is  the  same  as  the  thickness  6^  of 
the  plume  at  the  ceiling.  The  more  general  assumption  is  that  the  ceiling  jet  will  either  start  at 
or  expand  to,  some  width  Wj,  such  that  its  thickness  is  A (obtained  experimentally  or  via  some 
other  argument),  rather  than  6^.  We  then  have,  for  the  ceiling  jet. 


and  Eq.(4)  implies 


Fm^  = 


wv 

j c 


Pgc^ 


(13) 


Frh, 


Wj  = 


^cPg,c 


— = FwB.IA 
A ^ 


(14) 


We  therefore  need  to  find  F as  a function  of  w,  h^,  and  the  total  power  output  of  the  fire,  Q. 
We  can  express  it  in  terms  of  w by  using  Eq.(5);  then  Eq.(14)  becomes 


w-  = w — 
^ A 


1 

w 


= _f(>v-25) 
A 


(15) 


It  is  possible,  as  discussed  earlier,  that  the  flow  for  the  central  ceiling  jet  fans  out  radially  from 
its  initial  width  (at  the  ceiling)  of  Wj. 


In  order  to  discuss  the  phenomena  involved,  precisely,  refer  to  figure  14a,  a schematic  of  the 
appearance  of  wall-fire  flames  which  strike  Ae  ceiling  once  the  pyrolysis  zone  has  extended 
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beyond  z = h^.  The  ceiling  is  allowed,  in  imagination,  to  be  swung  upward  as  if  it  were  hinged 
at  its  intersection  with  the  wall.  The  pyrolyzing  areas  are  indicated  by  the  unprimed  letters,  and 
the  preheating  areas,  by  primed  letters.  Thus,  A,  B,  and  C are  the  pyrolyzing  areas  (and  A" 
as  well),  while  A’,  B’,  and  C’  are  preheating  areas.  Figure  15  shows  a few  more  details. 

Then  consider  the  case  illustrated  in  figure  4;  that  is, 

Xp  < h^  < Xf.  (16) 

Then  the  flames  in  Section  C’  (C  being  nonexistent  for  an  inert  ceiling)  irradiate  A",  both  in  the 
pyrolyzing  and  in  the  preheating  sections.  Thus,  enhanced  power  outputs  and  enhanced  spread 
rates  are  to  be  expected. 

Although  there  has  been  preheating  in  the  region  under  the  extensions,  that  heating  is  sufficient 
to  permit  lateral  growth  along  the  arms  (sections  B)  only  when  Xp  reaches  h^. 

We  next  need  to  discuss  the  view  factors.  This  is  sufficiently  complex  that  a separate  section 
is  devoted  to  it. 


4.3  Radiation  Feedback 

4.3.1  Experimental  data 

The  most  nearly  relevant  experimental  data  are  from  Hinkley  et  al  (1984).  Some  of  their  results 
are  given  in  Appendix  A.  Not  all  of  the  data  given  there  will  be  used  in  this  paper.  However, 
it  is  important  that  the  functional  form  of  the  fluxes  be  shown.  Moreover,  it  is  important  that 
we  note  the  observations  made  by  Hinkley  et  al  on  ignition  of  and  spread  of  fire  over  the  floor, 
as  an  indication  of  the  importance  of  the  ceiling  radiation  feedback.  Finally,  these  data  will  be 
used  in  our  subsequent  studies. 

With  their  configuration,  the  peak  flux  to  the  ceiling  was  170  kW/m^.  The  peak  downward  flux 
was  150  kW/m^;  this  is  high  enough  to  ignite  most  combustible  materials  within  a few  seconds. 
Hinkley  et  al  assumed  that  there  was  a wood  floor;  they  found  that  the  flux  was  sufficiently 
intense  that  the  floor  would  ignite  and  that  the  resulting  fire  would  spread.  They  carried  out  the 
calculation  for  several  assumed  ceiling  types,  including  a non-combustible,  perfectly  insulating 
ceiling;  an  asbestos  ceiling;  and  several  types  of  flammable  ceilings.  For  the  non-flammable  as 
well  as  for  several  of  the  flammable  ceilings,  the  results  were  very  similar  to  the  results  that  are 
observed  for  the  spread  of  flame  on  a vertical  wall:  for  wood  materials,  the  upward  spread 
ceases  beyond  a certain  point;  the  pyrolysis  front  attains  an  asymptotic  value.  (For  some  other 
flammable  ceilings,  the  downward  flux  is  so  high  that  spread  continues  and  even  accelerates). 
Whether  either  of  these  behaviors  could  occur  for  a vertical  flame  alone  is  doubtful. 
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4.3.2  Theoretical  calculations  of  fluxes 


The  radiation  feedback  from  the  arms  to  the  wall  (region  A"  in  figure  14b)  is  small,  since  the 
fiames  are  largely  in  the  same  plane.  Moreover,  their  contribution  is  not  very  different  from 
what  it  would  be  from  the  flames  above  A",  if  the  flames  extended  vertically,  as  in  figure  14b. 
Therefore  the  additional  fluxes  seen  by  the  pyrolyzing  and  preheating  wall  in  the  original  v 
ertical  section,  as  well  as  in  the  arms,  arise  mostly  from  the  flame  region  C,C’  on  the  ceiling 
(see  figure  14a).  The  radiation  from  A",  B,  and  B’  to  C and  C’  is  only  relevant  for  flammable 
ceilings.  Since  we  are  limiting  ourselves  to  inert  ceilings  for  the  present,  it  is  not  necessary  to 
calculate  them. 


The  radiation  feedback  from  the  flame  on  the  ceiling  (the  areas  C and  C’  in  Figure  14a,  or  the 
area  A2  in  figure  17) ) to  the  original  vertical  wall  sections  (A"  in  Figure  14b)  and  to  the  arms 
(areas  B,  B’  in  Figure  14a)  are  calculated  in  Appendix  B.  Note  that  with  the  assumption  that 
the  shaded  parts  of  A and  A’  in  figure  14b  correspond  to  B and  B’  in  figure  14a,  the  unshaded 
parts  of  A and  A’  correspond  to  C and  C’  in  14a.  The  power  supplied  to  the  lower  part  of  the 
wall  is  given  by 


TT 


tan-'i-Itan-'i.Iln 

(1  *r^)Y^ 

> 

1 

to 

7^(1  + 

(B25) 


where  corresponds  to  the  normalized  thickness  of  the  hot  layer  and  corresponds  to  the 
room  height.  These  and  the  other  symbols  are  defined  in  the  Appendix.  The  mean  radiation 
flux  from  the  ceiling-jet  flame  is  given  by 


4>,(0) 


(B38) 


where  <t>^{0)  is  the  downward  radiation  flux  from  the  flame  at  the  intersection  of  the  wall  and 
ceiling. 


The  power  supplied  to  the  upper  part  of  the  wall  is  given  by 


T 


, -1  1 r'  -1  1 y'  . 

tan  — - — tan  — + — Ini 

(1  +r'^)y'^ 

» 

y / j./  J.I  2 

(B30) 


where  Y’  = Y-J2  and  r’  is  given  by  Eq.(B32). 

Finally,  the  mean  flux  falling  at  the  distance  z downstream  on  each  of  the  arms  (see  figure  18) 
is,  from  Eq.(B36), 
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c 


(17) 


4>. 


arm 


^ ^rm 

tan  ^ 

z+b 

-tan  ^ 

z 

2t 

c 

c 

r •> 

tan"^ 

r 

1 

^ 

- tan  ^ 

1 

r 

'■i 

where  r is  given  by  Eq.(B22)  and  c and  z are  given  by  Eq.(B35). 


5.  THE  ALGORITHM 

The  most  attractive  approach  for  establishing  an  algorithm  for  the  growth  and  spread  of  the  fire 
along  wall  and  ceiling  would  be  to  use  a simple  set  of  algebraic  equations.  Unless  we  are 
merely  to  use  correlations,  however,  it  is  apparent  that  we  must  solve  coupled  differential 
equations  which  describe  various  rates.  Indeed,  a number  of  such  analytic  formulations  of  the 
upward  spread  process  have  been  developed  (see  Saito  et  al  1985;  Mitler  1986;  Quintiere  et  al 
1986).  The  most  recent  is  due  to  Cleary  and  Quintiere  (1991);  however,  they  generally  require 
simplistic  approximations  and  limitations  to  be  made.  The  alternative  is  to  use  a numerical 
approach,  which  can  be  more  general  and  more  powerful.  This  was  done  for  the  upward  spread 
process,  and  the  result  was  the  program  SPREAD  (Mitler,  1990),  which  is  described  in 
Appendix  C.  As  discussed  above,  one  of  the  limitations  of  SPREAD  is  that  it  does  not  take  the 
ceiling  into  account  except  in  a limited  way.  We  now  describe  how  to  overcome  this  deficiency. 


5.1  What  is  Needed  for  the  Algorithm 

As  indicated  in  the  discussions  above,  we  need  to  have  the  width,  thickness,  and  length 
(extension)  of  the  arms,  as  functions  of  the  several  independent  parameters.  The  same  is 
required  for  the  central  ceiling  jet  and  the  length  and  shape  of  that  part  of  the  flame.  The  heat 
fluxes  and  radiation  view  factors  from  the  parts  of  the  flames  lying  against  the  ceiling  back  to 
the  ceiling  and  wall  are  also  needed.  A procedure  was  previously  described  to  calculate  the 
various  view  factors.  Section  4.3.  Other  variables,  such  as  5,  y^,  and  a (in  Figures  15 

to  18)  must  be  determined  by  making  some  simple  assumptions. 


5.2  A First  Approximation 

Consider  the  growth  of  the  pyrolysis  region.  We  wish  to  transform  the  calculations  made  by 
SPREAD  for  flame  and  pyrolysis  front  positions  for  z > h^  into  the  actual  distributions  (over 
the  central  jet,  if  any,  plus  along  the  arms). 

Recall  that  Figure  14a  is  a schematic  of  the  appearance  of  the  fire,  for  the  case  that  the  pyrolysis 
zone  has  extended  beyond  z = h^.  The  areas  marked  B and  C are  the  pyrolyzing  areas,  while 
B’  and  C’  are  the  "preheating"  regions  where  flame  covers  the  surface  but  the  latter  has  not  yet 
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reached  the  pyrolysis  temperature.  For  a non-flammable  ceiling,  C has  zero  area,  so  that  the 
dashed  boundary  between  C and  C’  would  lie  precisely  at  the  wall-ceiling  intersection. 

Figure  14b  shows  the  geometry  implicit  in  SPREAD;  again,  A is  the  pyrolyzing  area  on  the 
ceiling.  A’  the  preheating  area,  and  A"  is  the  pyrolyzing  area  on  the  wall. 

For  the  zeroth-order  approximation,  we  assume  that  the  thickness  of  the  plume  flows  along  B,B’ 
and  C,C’  in  Figure  14a  are  exactly  the  same  as  for  A, A’  in  14b;  that  is,  in  the  notation  of 
Section  4.2,  6^  — A = Then  from  mass  conservation,  the  flow  divides  into  B and  C in 
proportion  to  their  respective  widths,  as  implied  by  Eq.(8).  In  other  words,  the  extensions  B,B’ 
are  just  the  same  as  the  shaded  parts  of  the  areas  A and  A’  in  flgure  14b,  of  width  5,  while  the 
region  C,C’  has  the  width  w - 25. 

We  have  assumed  that  the  wall  is  flammable  everywhere;  the  resulting  rapid  lateral  spread  of 
the  arms  across  the  wall  translates,  in  this  approximation,  into  continued  vertical  growth  beyond 
z = h^.  The  horizontal  extensions  (the  arms)  then  are  simply  the  shaded  sections  of  AH- A’  in 
Figure  14b,  turned  through  90°.  Thus  their  length  is 

= Xf  - Kv, 

and 

ye  = w/2  + Xp  - Iv 

(see  flgure  15,  which  complements  Figures  7 and  14). 

If  the  ceiling  is  flammable,  this  is  all  straightforward.  On  the  other  hand,  if  the  ceiling  is  inert, 
the  situation  is  apparently  inconsistent,  since  the  pyrolyzing  section  C has  zero  extent,  whereas 
B is  flnite.  That  is,  the  region  beyond  z = h^  bums  in  the  horizontal  orientation,  whereas  it 
cannot,  in  the  vertical  orientation  (since  the  ceiling  is  not  flammable):  the  shaded  parts  of  A and 
A’  in  flgure  14b  bum,  the  unshaded  parts  do  not.  The  procedure  which  resolves  this  is  to 
continue  to  calculate  the  flame  height  Xf  as  if  the  entire  region  were  in  fact  burning  ~ that  is, 
use  Q*  according  to 


(18) 

(19) 


y 

Q*  = wH^  j m\z)dz 


to  And  Xf,  but  then  note  that  the  actual  power  output  is  given  by 


y 

Q =wH^  f m"{z)dz  + f m"{z)dz 


(20) 


(21) 


That  is,  the  length  of  C is  zero,  that  of  C’  is  approximately  constant  at  Xf(t=tc)  - where  t^ 
is  the  moment  when  the  pyrolysis  front  on  the  wall  arrived  at  the  ceiling.  Thus  for  a non- 
flammable ceiling,  the  calculated  area  A must  be  multiplied  by  25/w  before  adding  that  pyrolysis 
contribution  to  the  rate  of  heat  release. 
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The  procedure,  then,  is  the  following:  we  specify  the  slab  height  to  be  substantially  greater  than 
h^.  We  then  calculate  Xp(t)  with  SPREAD  in  the  usual  way.  The  program  is  modified, 
however,  so  that  when  Xf  > h^,  we  begin  to  calculate  ^^(t),  according  to  Eq.(18)  (as  better 
expressions  are  found,  they  will  be  used  in  lieu.  We  must  also  calculate  6(t),  according  to 
Eqs.(7)  or  (11)  (or  another  expression  to  be  dictated  by  experiment,  in  the  future).  Finally, 
when  Xp  grows  beyond  h^,  we  must  modify  the  expression  for  Q(t),  the  time-dependent  rate  of 
heat  release.  This  quantity  is  normally  found  in  SPREAD  by  summing  the  mass-loss  rates  for 
all  the  horizontal  strips  into  which  the  slab  is  divided,  then  multiplying  the  sum  by  the  effective 
heat  of  combustion,  Xa^c-  For  Xp  > h^,  we  must  then  multiply  the  mass-loss  rate  riij  for  the 
strips  j corresponding  to  x > h^,  by  26/w,  then  sum. 

Now  consider  the  modifications  that  must  be  made  in  the  zeroth-order  process,  due  to  the  extra 
heating  fluxes  from  the  hot  layer,  hot  ceiling,  and  flame  feedback:  the  pyrolysis  rate  in  region 
A"  is  enhanced  by  the  hot  layer  fluxes  and  by  the  radiation  from  the  orthogonal  flame  covering 
C,  C’.  Likewise,  the  pyrolysis  rate  in  the  shaded  part  of  the  model  region  A is  enhanced  by 
radiation  from  C,C’  to  B,B’.  These  fluxes  are  all  relatively  easily  calculated  from  the 
expressions  found  above.  The  flames  covering  the  regions  B and  B’  contribute  just  about  as 
much  to  A"  as  the  vertical  parts  A, A’  would,  were  the  wall  indeed  very  high. 

On  the  other  hand,  consider  the  differences  in  the  two  cases,  for  a moment.  In  the  no-ceiling 
case,  buoyancy  continues  to  accelerate  the  gases  (up  to  a point),  whereas  it  d^elerates  in  the 
arms  (due  to  wall  friction).  For  this  and  other  reasons,  entrainment  into  the  arms  will  be  at  a 
different  rate.  Again,  the  heating  fluxes  from  those  flames  back  to  the  wall  will  probably  be 
somewhat  different  from  those  calculated  for  the  upward  spread  case.  Finally,  these  will  have 
to  be  consistent  with  experimental  data.  These  and  other  corrections  to  the  crude  approximations 
made  above,  will  be  made  next  year. 

If  there  exists  a critical  flow,  then  we  expect  something  like  Eq.(ll)  for  5,  whereas  if  there  is 
no  critical  flow,  Eq.(7)  would  be  more  appropriate.  When  there  will  be  relevant  experimental 
data  to  appeal  to,  better  expressions  for  f ^ and  for  b will  be  used. 

Thus,  the  procedure  is  to  increase  the  "external"  fluxes  impinging  on  A,  A’,  and  A" 
appropriately,  calculate  the  resulting  spread  rate  via  SPREAD,  take  the  appropriate  fractions  for 
the  arms,  and  we  are  (in  this  simplest  of  approximations)  done. 


6.  SUMMARY  AND  DISCUSSION 

It  has  been  demonstrated  that  the  effects  that  a ceiUng  has  on  the  progress  of  a room  fire  can  be 
quite  significant.  We  have  described  a number  of  these  effects.  It  has  been  demonstrated  that 
the  heating  effects  of  the  hot  layer  (and  hot  ceiling)  can  already  be  taken  into  account  with 
current  models. 

A "zeroth-order"  algorithm  has  been  formulated,  which  involves  simple  assumptions  about  the 
fluid  flow  and  some  straightforward,  albeit  complicated,  expressions  for  the  radiation  feedbacks. 
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It  has  been  found  that  this  can  be  implemented  with  an  extension  of  the  numerical  approach 
currently  used  in  SPREAD  in  order  to  describe  the  growth  and  spread  of  the  fire  under  the 
ceiling  is  straightforward.  We  have  also  looked  into  the  possibility  of  using  some  simple 
analytic  formulation  to  do  this,  but  such  a formulation  presents  a number  of  problems. 

Two  effects  have  not  been  considered  in  developing  this  algorithm:  combustion  and  spread  in 
an  oxygen-vitiated  environment,  such  as  is  found  under  the  ceiling,  and  the  effect  found  by  Saito 
(see  Section  3.4).  These  should  be  investigated. 
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NOMENCLATURE 


b width  of  radiating  (flame)  zone;  essentially,  w 
B pyrolyzing  area  in  the  arms 
d depth  of  the  layer  of  hot  gases 
D distance  of  the  burner  beneath  the  ceiling 
E■^  power  impinging  on  area  Aj 

F fraction  of  plume  mass  flow  at  the  ceiling  which  continues  as  a ceiling  jet 
Fj.j  view  factor  of  area  j from  area  i 
g acceleration  of  gravity 
h^  height  of  ceiling  (i.e.,  of  wall) 

H(x)  (total)  heat  flux  to  ceiling 
I(x)  downward  radiation  flux 

f flame  length  measured  from  the  virtual  origin,  in  Hinkley’s  experiments  (see  Eq.(A4)) 
fg  length  of  extension  (i.e.,  of  arm);  see  figures  9 and  14 
m’  mass  rate  of  flow  of  fuel  per  unit  width  of  corridor 
m mass-loss  rate;  mass-flow  rate 
m^  mass  flow  in  the  arms 

(upward)  mass  flow  in  the  plume,  at  the  ceiling 
p parameter  in  Eqs.(7)  and  (B9) 

Q power  output  of  flame  (RHR) 
r dimensionless  distance;  see  Eq.(B2) 
u reciprocal  of  r;  see  Eq.(B18) 

V upward  velocity  of  thermal  (wall)  plume 
Vj,  upward  velocity  of  plume  at  the  ceiling 
w width  of  wall  fire/plume 
Wj  width  of  ceiling  jet;  see  Eq.(14) 
w^  critical  width  of  vertical  plume 

X coordinate  indicating  vertical  position  (horizontal  position,  for  Eqs.(Al)  and  (A2) 

Xf  flame-tip  height 
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Xq  distance  of  a "virtual  origin"  behind  the  rear  of  a corridor  [Eqs.(Al)-(A3)] 
Xp  pyrolysis  front  height 
X dimensionless  distance;  see  Eq.(B2) 
y coordinate  for  horizontal  position  on  the  wall 
y^  extent  of  pyrolysis  zone  in  arm  (see  Figure  14) 

Y dimensionless  distance;  see  Eq.(B2) 

Yj  dimensionless  thickness  of  hot  layer 
Y^  dimensionless  height  of  room 

z coordinate  for  horizontal  position  orthogonal  to  the  wall 


Greek  symbols: 

6 width  of  arms  (extensions) 

A thickness  of  ceiling  jet;  see  Eq.(13) 

<j)j.  flux  from  ceiling  flame 

<f>^  mean  flux  from  ceiling  flame 

Pa  density  of  the  ambient  air 

Pg  density  of  plume  gas 

Pgj.  density  of  plume  gas  at  the  ceiling 

8 thickness  of  the  plume 

6^  thickness  of  arm,  where  it  meets  the  wall  plume 
6^  thickness  of  wall  plume  at  the  ceiling 
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Appendix  A.  Experimental  Data 


The  most  nearly  relevant  data  here  are  from  Hinkley  et  al  (1984).  They  had  a burner  against 
an  inert  wall,  with  a "corridor"  ceiling  against  which  the  flames  splashed.  The  downward 
radiation  flux  (from  hot  gas/flame  and  hot  ceiling)  was  found  to  decrease  exponentially  with 
distance  from  the  rear  of  the  corridor: 

/(x)  = 7^exp[-4.0(x +x^)//]  (Al) 

where  x is  the  distance  measured  from  the  rear  of  the  corridor  (in  meters),  x^  was  the  distance 
of  a "virtual  origin"  behind  the  rear  of  the  corridor,  and  i was  the  flame  length  measured  from 
the  virtual  origin.  is  the  effective  intensity  of  radiation  at  the  virtual  origin,  and  had  the  value 
of  500  kW/m^. 

The  flux  to  the  ceiling  was  also  measured;  it  consisted  of  the  convective  heating  flux,  as  well 
as  the  radiation  from  the  hot  gases/flame,  and  had  a similar  spatial  dependence: 

H{x)  = H^exp[-4.6(x+xJ//]  (A2) 

with  Hq  = 700  kW/m^.  With  their  conflguration,  the  peak  flux  to  the  ceiling  was  170  kW/m^. 
The  peak  downward  flux  was  150  kW/m^;  this  is  high  enough  to  ignite  most  combustible 
materials  within  a few  seconds.  Hinkley  et  al  assumed  that  there  was  a wood  floor;  they  found 
that  the  flux  was  sufficiently  intense  that  the  floor  would  ignite  and  that  the  resulting  fire 
would  spread.  They  carried  out  the  calculation  for  several  assumed  ceiling  types,  including  an 
asbestos  ceiling,  a non-combustible,  perfectly  insulating  ceiling,  and  several  types  of  flammable 
ceiling.  For  the  non-flammable  as  well  as  for  several  of  the  flammable  ceilings,  the  results  were 
very  similar  to  the  results  that  are  observed  for  the  spread  of  flame  on  a vertical  wall:  for  wood 
materials,  the  upward  spread  ceases  beyond  a certain  point;  the  pyrolysis  front  attains  an 
asymptotic  value.  (For  some  other  flammable  ceilings,  the  downward  flux  is  so  high  that  spread 
continues  and  even  accelerates).  Whether  either  of  these  behaviors  could  occur  for  a vertical 
flame  alone  is  doubtful. 

The  location  of  the  virtual  origin  correlated  very  well  with  D,  the  distance  of  the  burner  beneath 
the  ceiling: 

* 0.66  + 0.682D  + 0.641)2  (A3) 


The  flame  lengths  f (measured  from  the  virtual  origin)  were  considerably  greater  than  the 
analogous  flame  heights.  For  their  experiments,  using  "town  gas,"  they  found  that 


I = 220 


m 
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(A4) 
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Here  m’  is  the  mass  rate  of  flow  of  fuel  per  unit  width  of  corridor,  g is  the  acceleration  of 
gravity,  and  is  the  density  of  the  ambient  air.  This  expression  is  valid  for  f /d  < 24,  where 
d is  the  depth  of  the  layer  of  hot  gases. 

The  dependence  shown  in  Eq.(A4)  is  given  because  it  turns  out  to  be  precisely  the  usual 
correlation  found  for  wall  flres,  expressed  in  terms  of  mass  flow  rather  than  enthalpy  flow. 
Moreover,  the  coefficient  is  about  four  times  larger  than  in  the  latter  case;  that  is,  the  horizontal 
flames  were  much  longer  than  ones  rising  vertically  from  the  same  fuel  bed.  This  has  important 
implications  about  entrainment  and  combustion  which  will  be  pursued  later;  it  suffices  to 
emphasize  them,  at  the  moment. 

The  exponential  decrease  of  the  heat  fluxes  found  in  Eqs.(Al)  and  (A2)  is  interesting,  in  light 
of  the  fact  that  in  the  intermittent  region,  vertical  flames  have  a power-law  decrease  in  total  heat 
flux  (Mitler,  1990).  Kulkami  et  al  (1991,  1992),  however,  have  made  an  equally  good  fit  to 
those  fluxes  with  an  exponential  distribution,  and  this  is  completely  analogous. 

Again,  the  difference  between  (Al)  and  (A2)  arises  principally  because  the  latter  includes 
convective  heat  transfer.  Similar  exponentially-decreasing  fluxes  (radiative  and  total)  can  be 
expected  for  the  flame  extensions  being  considered  here. 
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Appendix  B.  View  Factors  and  Fluxes 

Consider  the  radiation  back  to  the  original  vertical  wall  section,  first  (see  A",  figure  13b,  and 
figure  15).  Thus,  the  radiation  from  the  fiame  on  the  ceiling  to  a strip  running  through  point 
P,  on  the  wall,  is  wanted.  The  radiation  fiux  from  the  flame  to  the  wall  at  any  height  is  given 
in  the  program  SPREAD;  it  is  assumed  that  the  outward  radiation  is  the  same  as  that  toward  the 
wall,  at  every  height.  For  z > h^,  we  now  have  that  outward  flux  radiating  back  to  the  floor, 
the  wall,  and  the  arms  (on  the  wall). 


Consider  the  area  A2,  of  length  a,  width  b,  radiating  to  the  narrow  strip  of  area  dA^,  located 
at  the  distance  c (see  figure  15).  We  have  available  the  view  factor  from  dA^  to  A2;  it  is  (Siegel 
and  Howell,  1981) 


where 


1 

tan  — + — log 

Y^(r^  + 1) 

y.  -1 1 

- — tan  _ 

TT 

Y 2 

r^(Y^  + l) 

r r 

X S a/b,  Y S c/b,  and  + Y^. 


(Bl) 

(B2) 


If  the  outward  flux  from  the  area  A2  were  uniform,  we  would  obtain  the  energy  flux  impinging 
on  dA^  by  using  the  view  factor  F^2-dAi  which  is  obtained  from  F^i_^2  the  reciprocity 
relation 


(B3) 


However,  the  emitted  radiation  varies  from  point  to  point,  and  the  power  impinging  on  dA^  is 
therefore  given  by 


dE^  - f dA2Fj2-di(^^ 
^0 


(B4) 


where  the  subscripts  dl  and  62  are  shorthand  for  dA^  and  dA2,  respectively. 


We  therefore  first  find  the  (differential)  view  factor  E^i_^2  ^ strip  of  length  b and  width  da, 

at  a,  radiating  to  the  strip  (of  width  dc)  at  P. 


Since  dA2  = b da  (where  a=bX)  and  dAj  = b dc,  Eq.(Bl)  implies  that 

,-ii 


F,  , =dx^-  ax^Ym-\Vr) 


(B5) 


From  Eq.(B4)  and  the  reciprocity  theorem,  we  have 
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dE^  - I 


(B6) 


and  using  equation  (B5),  (B6)  becomes 

+ (B7) 

■'o  {X^  + y2)3/2 

This  integral  can  be  done  exactly  for  the  case  <}>j.  = const.  On  the  other  hand,  we  would  expect 
that  the  flux  falls  off  exponentially,  as  given  by  Eq.(18).  For  an  exponential  or  a power-law 
fall-off,  however,  the  integral  in  Eq.(28)  is  not  analytically  tractable.  But  if  we  use  the  mean- 
value  theorem  to  write  the  integral  in  Eq.(28)  as 


M = tan'i  {X^  + 


then  with  the  choice 


<t>,W  = <l),(0)e-«"  = 4>,(0)€-«»*, 


(B8) 


(B9) 


for  as  given  by  Eq.(18),  an  exact  integration  can  be  carried  out.  This  results  in 

(BIO) 

where  Nq  is  the  zeroth-order  Neumann  function  (identical  to  Yq,  the  Bessel  function  of  the 
second  kind),  Hq  is  the  Struve  function  of  order  zero,  and 

11  s qbY 


dEi  = 


^[tan-i(x2  + |l  ^ l7[iVo(7) -//(,(,,)] 


A calculation  shows  that  the  mean  value  of  X is  about  Y/2,  and  that  yields 

0.9 


. dA-i  _ . 

dE.  = tan  ^ 

TT 


<#.,(0)  Jl  + lv[No{v)-Ho(v)] 


(B12) 


We  can  obtain  a result  in  terms  of  more  familiar  functions  with  the  choice 


for  (this  p is  unrelated  to  the  one  found  in  Eq.(7));  that  yields 


dA-l  . — fy 

— L[tan"^(X^ 


where 


(B14) 
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S ^ \fpbY 


(B15) 


A simple  alternative  is  to  assume  that  the  radiation  flux  from  A2  has  some  mean  value  out 
to  the  flame-tip  height  and  zero  beyond.  In  that  case,  with  the  deflnitions 


and 


a = Xf  - hw, 

= a^  = (Xf  - h„)/b, 


u = 


1 


b 


+ Xl 


2 2 
c + a 


(BI6) 

(B17) 

(B18) 


Eq.(B7)  yields  Eq.(B19): 
dA 


dE^  = 


IT 


Y<t> 


rm 


itan  ^(1/y)  - Mtan  + 0.5  log(l  + + 0.5  log 


1 + 

While  this  eliminates  the  need  to  make  a judicious  choice  for  the  parameter  p in  Eq.(Bll),  it 
requires  an  equally  judicious  choice  for  the  mean  value  of  0^,  in  Eq.(B19).  On  the  other  hand, 
this  approach  also  permits  the  radiation  to  the  arms  to  be  readily  found,  as  well.  Note  that 
Eq.(B19)  is  just  Eq.(Bl). 

Next,  we  want  the  radiation  impinging  on  the  arms,  from  the  ceiling  flame.  The  power  to  the 
arms  is  again  given  by  Eq.(B6),  but  with  the  appropriate  view  factor.  The  view  factor  from  a 
surface  element  dAj  on  the  wall,  a distance  c down  from  the  ceiling,  to  the  flnite  rectangle  A2 
on  the  ceiling  (see  Figure  16)  is 
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with  the  same  definitions  for  X,  Y,  and  r as  in  Eq.(B2).  Then  by  a simple  subtraction,  we  find 
that  from  dA^  to  the  shaded  area  A3  in  figure  17,  we  have 
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where 
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and 


r 
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(B23) 


Note  that  we  can  non-dimensionalize  z to  Z by  dividing  z by  b,  as  was  done  for  X and  Y;  this 
is  not  necessary  here,  however.  If  we  again  assume  a uniform  flux  to  be  emitted  from  A3,  then 
the  energy  falling  on  the  area  dA^  on  the  arm  is 

dE,  = , (B24) 


where  we  have  again  used  the  reciprocity  theorem. 

Thus,  we  can  get  the  needed  fluxes  to  the  wall  and  to  the  arms  from  Eqs.(B19)  and  (B24), 
respectively.  At  present,  however,  SPREAD  assumes  a uniform  external  flux  on  the  section  of 
wall  within  (i.e.,  contiguous  to)  the  hot  layer  and  another  to  the  layer  below  ("upper"  and 
"lower  wall",  respectively).  Therefore  Eq.(B19)  must  be  integrated  between  z = z^  (the 
interface  height)  and  z = h^  (the  ceiling  height),  and  Eq.(B24)  between  h^  - 5 and  h^. 


Integration  of  Eq.(B19)  over  Y is  tedious  but  straightforward,  and  yields 
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and 
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We  can  again  simplify  this  very  considerably  by  taking  an  average  value  for  Y in  Eq.(B19). 
Thus,  starting  again  from  Eq.(B19)  we  find,  for  the  lower  wall. 
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where 

Yi  = ih^-Zi)/b 

corresponds  to  the  normalized  thickness  of  the  hot  layer,  and 


Y„.  = h,/b 


corresponds  to  the  room  height.  Also, 
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We  next  turn  to  Eq.(B24).  We  again  let  dA^  represent  the  entire  area  of  interest.  Thus  for  the 
pyrolyzing  area  B (see  figures  13  and  14),  dAj  becomes 

B = (ye  - w/2)5,  (B33) 


while  for  the  preheating  region  B’,  dA^  becomes 

B’  = (f  e - Ye  + w/2)5.  (B34) 

Moreover,  choosing  the  midpoints  of  these  regions  to  be  representative,  we  get 

c = 6/2  and  z = [ye  ■ w/2]/2  (B35) 

so  that  Eq.(B6)  yields 

El  ^ ^^36) 

with  F^j^.3  given  by  Eq.(B21),  with  c and  z given  by  Eq.(B35). 

The  simplest  approximation  to  is  found  by  taking 
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Suppose  the  radiation  fiux  is  given  by  an  equation  such  as  (Al);  that  is,  by  Eq.(B9),  where  0/0) 
is  the  downward  radiation  fiux  from  the  flame  at  the  intersection  of  the  wall  and  ceiling.  Then 
the  integral  on  the  left-hand  side  of  Eq.(B37)  is  trivial,  and  we  find  that 
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Appendix  C.  The  program  SPREAD 


The  following  is  a brief  description  of  SPREAD:  For  a given  igniter  against  a flat  wall,  and 
given  the  appropriate  thermophysical  parameters  of  the  wall  material,  SPREAD  calculates  the 
time  to  ignition  of  the  wall,  the  location  and  extent  of  the  initial  pyrolyzing  area,  the  mass-loss 
rate,  the  rates  of  upward  and  lateral  spread  of  the  fire,  and  the  resulting  rate  of  heat  release. 

The  input  includes  the  experimental  rate  of  heat  release  from  the  Cone  Calorimeter  (this  is 
proportional  to  the  mass-loss  rate).  Thus,  the  calculation  automatically  takes  charring  and 
transient  pyrolysis  into  account.  The  input  also  includes  the  external  flux  impinging  on  the  wall 
as  a function  of  time  and  position,  which  the  program  then  takes  into  account  in  making  the 
above  calculations  (the  present  program  allows  for  just  tw  external  fluxes:  the  mean  external 
flux  impinging  on  the  part  of  the  wall  contiguous  to  the  upper  layer,  and  that  wetting  the  lower 
layer.  Generalization  to  more  than  two  fluxes  is  straightforward). 

For  the  lateral-spread  calculation,  the  program  takes  into  account  the  fact  that  the  upper  part  of 
the  wall  has  a different  temperature  than  the  lower  part,  so  that  the  lateral-spread  rate  is  greater 
in  the  upper  region.  Thus,  the  effect  of  having  two  layers  is  explicitly  taken  into  account. 

The  power  output  of  the  igniter  is  user-specified,  and  can  be  an  (almost)  arbitrary  function  of 
time.  Among  the  things  the  program  does  not  yet  have  in  it  is  the  effects  the  presence  of  the 
ceiling  has  on  the  burning  and  spread  rates  (except  insofar  as  the  ceiling  traps  the  upper  gas 
layer). 

The  performance  of  SPREAD  (before  the  changes  indicated  in  this  paper  were  added)  for  three 
examples  is  illustrated  in  figures  19  to  21.  Figure  19  shows  the  position  of  the  pyrolysis  front 
for  a PMMA  slab  ignited  at  the  bottom,  as  a function  of  time.  The  squares  show  the  measured 
positions  of  the  front,  the  solid  curve  the  value  calculated  by  SPREAD,  the  BFRL  model  which 
is  being  prepared  for  insertion  into  HAZARD.  Figure  20  does  the  same  for  a slab  of  wood. 
Figure  21a  shows  the  measured  values  of  pyrolysis  front  and  flame-tip  height  for  a large  slab 
of  plywood  preheated  for  20  min.  by  a 5 kW/m^  radiation  source,  then  ignited  at  the  bottom. 
Figure  21b  shows  the  corresponding  calculated  values. 

The  documentation  of  SPREAD  is  given  in  much  more  detail  in  Mitler  and  Steckler  (1994). 
Finally,  the  results  of  implementing  the  algorithm  are  shown  in  Fig.  22:  this  shows  the  results 
of  burning  a 0.3  m-wide  PMMA  slab  of  2.43  (S’)  height  in  the  open,  vs.  that  in  an  enclosure. 
The  rate  of  heat  release  is  shown.  Only  the  splashing  effect  and  the  spread  of  the  horizontal 
extensions  under  the  ceiling  are  included  --  the  effect  of  the  extra  flux  from  the  hot  layer  is  not 
shown.  It  was  assumed  that  the  mean  flux  from  the  flame  under  the  ceiling  is  5 kW/m^. 
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FIGURE  CAPTIONS 


Figure  1. 

Schematic  showing  how  the  hot  layer  trapped  under  the  ceiling  (as  well  as  the 
ceiling  itself)  provides  an  additional  external  flux  which  increases  the  mass-loss  rate  in 
the  pyrolyzing  region  and  increases  the  heating  rate  (and  therefore  the  spread  rate)  in  the 
preheating  region. 

Figure  2. 

Schematic  for  what  one  might  imagine  is  the  appearance  of  a wall  flame  so  high 
that  it  "splashes"  on  the  ceiling.  The  flgure  is  shown  as  if  the  ceiling  were  hinged  where 
it  meets  the  wall,  and  rotated  back  so  as  to  lie  in  the  same  plane  as  the  wall.  All 
analogous  figures  shown  hereafter  use  the  same  device. 

Figure  3. 

Schematic  of  what  is  actually  observed  (see  figure  5).  The  dotted  curve  indicates 
the  possible  shape  of  the  part  of  the  flame  which  lies  against  the  ceiling;  the  actual  shape 
is  actually  unknown:  it  is  generally  hidden  by  smoke. 

Figure  4. 

Schematic  showing  how  flame  extensions  in  the  direction  of  the  ceiling  jet  produce 
enhanced  radiation  feedback  to  the  wall. 

Figure  5. 

Geometry  assumed  by  Cleary  and  Quintiere,  analogous  to  figure  3.  The  shapes 
on  the  ceiling  are  semicircles.  Several  are  shown,  to  indicate  the  spread  of  the  fire. 

Figure  6. 

Photograph  of  a wall  fire  with  the  flame  impinging  on  the  ceiling.  From  an  experiment 
carried  out  at  CFR/NBS  by  B.  Lee  in  1986.  Photo  courtesy  of  R.  Peacock. 

Figure  7. 

Schematic  of  an  intermediate-width  fire,  indicating  the  gas  flows  (dashed  lines), 
the  flame  extensions,  and  some  of  the  relevant  dimensions. 

Figure  8. 

Calculated  RHR  from  a slab  of  PMMA  ignited  at  the  bottom,  as  a function  of 
time.  Curve  NC  corresponds  to  burning  in  the  open,  curve  C shows  the  result  of 
subjecting  the  slab  to  the  extra  heating  flux  from  the  hot  gases  trapped  in  an  ASTM 
room.  The  arrows  at  the  top  indicate  the  moments  at  which  the  pyrolysis  front  reaches 
the  top  of  the  slab,  in  the  two  cases. 

Figure  9. 

See  caption  for  Fig. 8.  The  flame  tips  are  computed  to  reach  the  ceiling  at  t = 
422  s.  The  dashed  curve  indicates  schematically  what  the  further  development  of  the 
RHR  might  look  like,  with  the  flames  shooting  out  along  the  arms.  The  growth  rate 
would  be  still  more  rapid  if  the  ceiling  were  flammable. 

Figure  10. 

For  a very  wide  fire,  the  ceiling  flow  is  nearly  two-dimensional. 

Figure  11. 

Schematic  for  a narrow  wall  fire:  flow  occurs  principally  along  the  arms.  The 
central  jet,  if  there  is  one,  may  proceed  as  a "tendril"  of  gas,  or  fan  out  radially,  as 
shown  by  the  dashed  curves. 

Figure  12. 

The  dependence  of  the  ceiling-jet  fraction  F,  as  a function  of  the  wall-fire  width. 
Curve  A:  for  A < l/(2p)  in  Eq.(8) 

Curve  B:  for  A = l/(2p) 

Curve  C:  when  there  exists  a critical  width,  w^. 
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Figure  13. 

The  width  8 of  the  arms  as  a function  of  the  width  w of  the  wall  fire,  according 
to  Eq.(l)  (dashed  curves)  and  Eq.(ll)  (solid  curve). 

Curve  A:  2Ap  < 1 

Curve  B:  2Ap  > 1 

Figure  14a. 

Schematic,  similar  to  figure  7.  C and  C’  are  the  luminous  (flame)  regions  within 
the  central  jet.  C is  the  pyrolyzing  region,  when  the  ceiling  is  flammable  (it  has  zero 
area  when  the  ceiling  is  inert).  B and  B’  are  the  analogous  regions  in  the  arms. 

Figure  14b. 

A and  A’  are  the  zones  for  the  virtual  wall  fire  which  are  used  to  model  the  actual 
fire;  that  is,  from  which  we  can  infer  the  fire  strength  and  spread  rate  for  configuration 
(a).  A"  is  the  actual  pyrolyzing  area  on  the  wall.  The  shaded  areas  are  explained  in  the 
text  (Section  5.2). 

Figure  15. 

This  figure  combines  figures  7 and  14a,  and  shows  a few  more  relevant 
dimensions. 

Figure  16. 

Schematic  showing  the  geometry  for  the  calculation  of  the  view  factor  at  point  P 
and  for  the  strip  running  through  P,  from  the  luminous  area  A2. 

Figure  17. 

Configuration  similar  to  that  shown  in  figure  16,  but  applying  to  the  differential 
element  dAj  at  the  edge  of  region  A2,  rather  than  the  entire  strip. 

Figure  18. 

Applying  the  view  factors  found  from  figure  17  twice,  we  find  the  view  factor  for 
a point  on  the  wall  under  the  arm,  located  a distance  z "downstream"  from  A2  and  a 
distance  c down  from  the  ceiling. 

Figure  19. 

Position  of  the  upper  pyrolysis  front  as  a function  of  time,  for  thermally-thick 
PMMA,  NIST  experiment.  The  difference  between  peak  and  average  values  arises 
because  the  pyrolysis  front  is  not  horizontal  (three-dimensional  burning).  The  solid  curve 
is  the  calculation. 

Figure  20. 

Position  of  the  upper  pyrolysis  front  as  a function  of  time,  for  particle  board. 
Same  remarks  as  for  Figure  19. 

Figure  21a. 

Summary  of  experiment  #17  (particle  board,  ~ ^ kW/m^). 

Figure  21b. 

Calculated  values  of  pyrolysis  front  height  (solid  curve)  and  of  flame-tip  heights 
(dotted  curve). 

Figure  22. 

Rate  of  heat  release  (in  kW)  of  a PMMA  slab  as  a function  of  time,  with  and  without  a 
ceiling. 
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Hot  Ceiling 


Figure  1.  Schematic  showing  how  the  hot  layer  trapped  under  the  ceiling  (as  well  as  the 
ceiling  itself)  provides  an  additional  external  flux  which  increases  the  mass-loss 
rate  in  the  pyrolyzing  region  and  increases  the  heating  rate  (and  therefore  the 
spread  rate)  in  the  preheating  region. 


32 


Ceiling,  Folded  Up 
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Figure  3.  Schematic  of  what  is  actually  observed  (see  figure  5).  The  dotted  curve  indicates 
the  possible  shape  of  the  part  of  the  flame  which  lies  against  the  ceiling;  the 
actual  shape  is  actually  unknown:  it  is  generally  hidden  by  smoke. 
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Figure  4.  Schematic  showing  how  flame  extensions  in  the  direction  of  the  ceilin, 
enhanced  radiation  feedback  to  the  wall. 
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Figure  5.  Geometry  assumed  by  Cleary  and  Quintiere,  analogous  to  figure  3 . The  shapes 
on  the  ceiling  are  semicircles.  Several  are  shown,  to  indicate  the  spread  of  the 
fire. 
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Figure  6. 


Photograph  of  a wall  fire  with  the  fiame  impinging  on  the  ceiling.  From  an 
experiment  carried  out  at  CFR/NBS  by  B.  Lee  in  1986.  Photo  courtesy  of  R. 
Peacock. 
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Figure  7.  Schematic  of  an  intermediate- width  fire,  indicating  the  gas  fiows  (dashed  lines), 
the  fiame  extensions,  and  some  of  the  relevant  dimensions. 
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with  ceiling  no  ceiling 
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Figure  8.  Calculated  RHR  from  a slab  of  PMMA  ignited  at  the  bottom,  as  a function  of 
time.  Curve  NC  corresponds  to  burning  in  the  open,  curve  C shows  the  result 
of  subjecting  the  slab  to  the  extra  heating  flux  from  the  hot  gases  trapped  in  an 
ASTM  room.  The  arrows  at  the  top  indicate  the  moments  at  which  the  pyrolysis 
front  reaches  the  top  of  the  slab,  in  the  two  cases. 
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Figure  9.  See  caption  for  Fig.  8.  The  flame  tips  are  computed  to  reach  the  ceiling  at  t 

422  s.  The  dashed  curve  indicates  schematically  what  the  further  development 
of  the  RHR  might  look  like,  with  the  flames  shooting  out  along  the  arms.  The 
growth  rate  would  be  still  more  rapid  if  the  ceiling  were  flammable. 


Figure  10.  For  a very  wide  fire,  the  ceiling  fiow  is  nearly  two-dimensional. 
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Ceiling 


Figure  11.  Schematic  for  a narrow  wall  fire:  flow  occurs  principally  along  the  arms.  The 
central  jet,  if  there  is  one,  may  proceed  as  a "tendrir'  of  gas,  or  fan  out  radially, 
as  shown  by  the  dashed  curves. 
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Figure  12.  The  dependence  of  the  ceiling-jet  fraction  F,  as  a function  of  the  wall-fire  width. 
Curve  A:  for  A < l/(2p)  in  1^.(8) 

Curve  B:  for  A = l/(2p) 

Curve  C:  when  there  exists  a critical  width,  w^ 


Figure  13.  The  width  6 of  the  arms  as  a function  of  the  width  w of  the  wall  fire,  according 
to  Eq.(l)  (dashed  curves)  and  Eq.(ll)  (solid  curve). 

Curve  A:  2Ap  < 1 
Curve  B:  2Ap  > 1 
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Figure  14a.  Schematic,  similar  to  figure  7.  C and  C’  are  the  luminous  (flame)  regions  within 
the  central  jet.  C is  the  pyrolyzing  region,  when  the  ceiling  is  flammable  (it  has 
zero  area  when  the  ceiling  is  inert).  B and  B’  are  the  analogous  regions  in  the 
arms. 


Figure  14b.  A and  A’  are  the  zones  for  the  virtual  wall  fire  which  are  used  to  model  the  actual 
fire;  that  is,  from  which  we  can  infer  the  fire  strength  and  spread  rate  for 
configuration  (a).  A"  is  the  actual  pyrolyzing  area  on  the  wall.  The  shaded  areas 
are  explained  in  the  text  (Section  5.2). 


44 


45 


Figure  15.  This  figure  combines  figures  7 and  14a,  and  shows  a few  more  relevant 
dimensions. 


Figure  16.  Schematic  showing  the  geometry  for  the  calculation  of  the  view  factor  at  point  P 
and  for  the  strip  running  through  P,  from  the  luminous  area  A2. 
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Figure  17. 


Figure  18. 


Configuration  similar  to  that  shown  in  figure  16,  but  applying  to  the  differential 
element  dAj  at  the  edge  of  region  A2,  rather  than  the  entire  strip. 


Applying  the  view  factors  found  from  figure  17  twice,  we  find  the  view  factor  for 
a point  on  the  wall  under  the  arm,  located  a distance  z "downstream”  from  A2 
and  a distance  c down  from  the  ceiling. 
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UPWARD  FLAMESPREAD  ON  PMMA 
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Figure  19.  Position  of  the  upper  pyrolysis  front  as  a function  of  time,  for  thermally-tliick  PMMA;  NIST 
experiment.  The  difference  between  peak  and  average  values  arises  because  tlie  pyrolysis 
front  is  not  horizontal  (three-dimensional  burning).  The  solid  curve  is  calculated. 


UPWARD  FLAMESPREAD 
CFR  Particle  Board 


CO 

CD  LO 

11  CN 

"x  J 
0 jCi 

a o 


o 

o 

o 


o 

o 


o 

o 

CD 


O 

O 

LO 


o 

o 


o 

o 

ro 


o 

o 

CNJ 


o 

o 


o 


o 

o 

d 


00 

0 


h- 


[UJ]  dx 


49 


Figure  20.  Position  of  the  upper  pyrolysis  front  as  a function  of  time,  for  particle  board. 
Same  remarks  as  for  Fig.  19. 


Figure  21b.  Calculated  values  of  pyrolysis  front  height  (solid  curve)  and  of  flame-tip  heights 
(dotted  curve). 
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Figure  22.  Rate  of  heat  release  (in  kW)  of  a PMMA  slab  as  a function  of  time,  with  and 
without  a ceiling. 
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